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Abstract 

An Ising spin system under the critical temperature driven by a dichotomous Markov noise (mag- 
netic field) with a finite correlation time is studied both numerically and theoretically. The order 
parameter exhibits a transition between two kinds of qualitatively different dynamics, symmetry- 
restoring and symmetry-breaking motions, as the noise intensity is changed. 

There exist regions called channels where the order parameter stays for a long time slightly above 
its critical noise intensity. Developing a phenomenological analysis of the dynamics, we investigate 
the distribution of the passage time through the channels and the power spectrum of the order 
parameter evolution. The results based on the phenomenological analysis turn out to be in quite 
good agreement with those of the numerical simulation. 
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I. INTRODUCTION 



Over the last decade, the dynamics of ferromagnetic systems below their critical tem- 
peratures in a periodically oscillating magnetic field have been studied both theoreti- 
cally and experimentally 10]. The systems exhibit two qualita- 
tively different behaviors referred to as symmetry-restoring oscillation (SRO) and symmetry- 
breaking oscillation (SBO), depending on the frequency Q and the amplitude h of the applied 
magnetic field. It has been established that there exists a sharp transition line between SRO 
and SBO on the (Q, h) plane, which is called the dynamical phase transition (DPT). The 
DPT was first observed numerically in the deterministic mean-field system for a ferromagnet 
in a periodically oscillating field l|, and has subsequently been studied in numerous Monte 
Carlo simulations of the kinetic Ising system below critical temperature 



also been observed experimentally in an ultra-thin Co film on Cu(100) 



It has 



Recently, we investigated the DPT by introducing the model equation s(t) 
T)s — s 3 + hcosQt 7]. This equation is a simplified model for the Ising spin system at 
the temperature T below its critical value T c in an external periodic magnetic field. By 
appropriately scaling the magnetization s, time t, and the applied field, this equation is 
written as 

s(t) = s - s 3 + hcosttt. (1) 

The SBO and SRO are observed in Eq. (JTJ) and the transition line between them on the (Q, 
h) plane is determined analytically ?], Isj]- 

It is quite interesting to ask whether DPT is observed under another kind of applied field, 
especially random field with bounded amplitude. The fundamental aim of the present paper 
is to study the dynamics of s(t) with a dichotomous Markov noise (DMN) F(t) instead of 
periodically oscillating external field hcosQt (see, e.g., [ll|). 

The equation of motion 

s = f(s)+F(t) (2) 
with a nonlinear function f(s) and the DMN F(t) has been extensively studied by many 



authors 



3,0, 0,0 



It is well known that the master equation for the system can 



be derived, and then transition phenomena of stationary probability densities concerning 
the intensity of F(t), for example, are studied, which are referred to as the noise-induced 
phase transition 12j, [l6[. The asymptotic drift velocity (s) in the case of f(s) being periodic 



functions are also discussed as a specific dynamic property Furthermore, the mean first- 
passage time (MFPT) and transition rates are investigated as another important dynamic 
property when f(s) is the force associated with the bistable potential given by Eq. (J2J) [3, 



151 ]. For a review of works on DMN system, see Bena [lj 

The fundamental aim of the present paper is to propose a phenomenological approach 
to the critical dynamics near the transition point between the symmetry-restoring motion 
(SRM) and the symmetry-breaking motion (SBM) observed in Eq. (J2J) with f(s) = s — s 3 
(Sec. |TTJ) and to investigate the distribution of passage times through channels, switching 
times between two different motions, and the power spectrum of the order parameter evo- 
lution. The present paper is constructed as follows. In Sec. |H1 we discuss the dynamics of 
symmetry-breaking motion and symmetry-restoring motion of the model equation (J2J). In 
Sec. IIII1 the jumping process of the magnetization through channels, which are defined in the 
text, is investigated and the MFPT is obtained. In Sec. IIVI a phenomenological approach 
simplifying the dynamics of passing through a channel in the SRM phase is introduced and 
three statistical characteristics are analytically developed. The results are compared with 
numerical simulations. Concluding remarks are given in Sec. 

II. MODEL EQUATION AND SYMMETRY-BREAKING TRANSITION 

A. Model equation and noise-induced phase transition 

We consider the equation of motion driven by the external field F(t), 

f(s) + F(t), (f(s) = s-s s ) (3) 



dt 

where F(t) is a symmetric DMN with taking the values ±Hq. Here the probability p{r) that 
Fit) continues to take the identical value +H or — H longer than time r is given by 

p{r)=e- T ' T f. (4) 

This implies that the correlation time of F(t) is equal to T//2. Throughout this paper, 
numerical integrations of Eq. (JHJ) are carried out by using the Euler difference scheme with 
the time increment At = 1/100. 

Without DMN, s(t) eventually approaches either of the stationary fixed points ±1, one of 
which is achieved according to the initial condition s(0) as shown in Fig. ^ In the presence 


























.r 




f 



















FIG. 1: The function f(s) = s — s 3 is shown by solid lines, (a) Real roots Sj±, (j = 1, 2, 3) for 
Ho < H c and (b) real roots s± for Hq > H c of the algebraic equation s — s 3 ±Hq = 0. The definition 
of H c is graphically represented. Two bold lines drawn in (b) indicate the channel regions defined 
in Sec. EH 



of DMN, if H < H c , H c being denned by 

H c = 2(l/3) 3/2 = 0.3849- • • 



(5) 



then f(s)+H = 0(f(s) — H Q = 0) has three real roots s J - + (s J -_), (j = 1, 2, and 3). Each value 
of Sj± is graphically shown in Fig. QJa). On the other hand, if H > H c , then f (s) + H = 
(f(s) — Ho = 0) has only one real root s + (s_) given by 



8± 



^0 2 - ^c 2 



nl/3 



+ 



±H n 



H 2 o ~ HI 



1/3 



(6) 



which are indicated in Fig. QJb). 

Next let us consider the dynamics described by Eq. (JSj) for H < H c and for H > H c , 
and discuss similarity and difference between the dynamics in the periodically oscillating 
field case and those in the present DMN case. A part of our results belon gs t o the context 
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In the case 



of the noise-induced phase transition and MFPT in Refs. 
of Ho < H c , three motions numerically integrated are shown in Figs. Efa) and (b). Two 
motions confined in the ranges si_ < s(t) < Si + and S3_ < s(t) < 53+ are both stable. The 
long time average (s(t)) of each motion does not vanish, and the motion is called SBM in 
relation to DPT in the oscillating external field case. On the other hand, the motion s u (t) 
confined in the range S2+ < s u (t) < S2- is unstable. The long time average of s u (t) vanishes, 
and in this sense the motion is called SRM. It should be noted that this unstable SRM is 
located between two stable SBM, which has a similar characteristic to SBO of DPT 3]. 
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FIG. 2: Figures (a) and (b) show the motions obtained by numerically integrating Eq. (J3J) for 
Hq = 0.38(< H c ) and tj = 5, where two SBM's (solid line) and an unstable SRM (dashed line) 
are drawn. The unstable SRM is evaluated by replacing t — > —t. On the other hand, Figs, (c) and 
(d) show the motions for Hq = 0.5(> H c ) and tj = 10. 

The motion of s(t) for Hq > H c is shown in Figs. Efc) and (d). One observes that there 
exists a stable SRM confined in the range s_ < s(t) < s + . For SRM, the time average of 
s(t) vanishes, i.e., (s(t)) = 0. The comparison between Figs. I2fb) and (d) suggests that the 
SRM for Hq > H c is generated via the "attractor merging crisis" [18| of the two SBM's and 
one unstable SRM, i.e., the two SBM's and one unstable SRM disappear and then one stable 
SRM takes place at H Q = H c . This situation is similar to that in the DPT case. However, 
in contrast to the DPT case, as will be shown in Sec. Ill B[ the transition line on the (tJ 1 , 
Hq) plane is independent of the correlation time r/ of F(t) and the average (s(t)) depends 
discontinuously on H . 



B. Stationary distribution functions and phase diagram 

In this subsection, we discuss the stationary distribution functions for SBM and SRM. To 
this aim, we first consider a slightly general nonlinear Langevin equation of motion driven 
by DMN, 

x(t) = f(x)+g(x)F(t), (7) 

where f(x) and g(x) are generally nonlinear functions of x and F(t) is DMN The 
temporal evolution of the distribution function P(x,F,t) that x(t) and F(t) respectively 
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take the values x and F(= ±H ) is determined by |l2t ll 

^P(x,t) = ~[f(x)P(x,t)+H g(x)q(x,t)], 
d 2d 

— q(x,t) = -—q(x,t)- — [f(x)q(x,t) + H g(x)P(x,t)}, (8) 

where we put P(x, t) = P(x, +H , t)+P(x, —H , t) and q(x, t) = P(x, +H , t)—P(x, —H , t). 
The stationary distribution P st (x) = P(x, oo) is solved to yield 



X 



7St( \ _ AT g ^ X > o XD j \ 

{) Hlg{xY-f{xY P \ T fJ ~~ [f(x')-H g(x') ' f(x>) + H g(x')\ 



1 1 

+ 



(9) 



provided that each of the equations 

x = f(x)+H g(x), x = f(x)-H g(x) (10) 

has at least one stable fixed point, where N is the normalization constant. 

By substituting f(x) = x — x 3 and g(x) = 1, (Eq. (J3J), into Eq. Q, the stationary 
distribution function P§bm( s ) f° r SBM (H < H c ) for s%- < s < s 3+ or si_ < s < si+ is 
written as 

PJW*) OC |« 2 - S ? + |-^|5 2 - .Ll^-k 2 - sl + \~^, (11) 

/5i± = l-T7 1 |(s J ±-Sfe±)(s i± -^±)|, (12) 

where (j, k, I) = (1, 2, 3), (2,3,1), and (3,1,2). On the other hand, the stationary distribution 
function Psrm( s ) f° r ^ ne SRM (if > H c ) for s_ < s < s + is obtained as 

-i -i 

-J 1 



-i. 



x exp 



— arctan 



(4-1/3)734^1 

2s + s+ 



arctan 



2s - s_ 



3sf 



(13) 



v/34^4 

The analytic solutions (jll|) and ([13]) are numerically confirmed in Fig. [HI 

As ifo is increased, the form of the stationary distribution function changes drastically 
from the forms in Eq. (fTT[) to Eq. (fTHj) at Hq = H c . This phenomenon which is induced 
by the disappearance of two pairs of stable and unstable fixed points [13] is an example 
of the noise-induced phase transitions 1(|. It turns out that the transition line between 




FIG. 3: Stationary distribution functions obtained both analytically (solid line) and numerically 
(symbols) for (a) Hq = 0.3 and tj = 5, and for (b) Hq = 0.39, 0.4, and 0.5 by keeping tj = 10. 
(a) and (b) correspond respectively to SBM and SRM. The distribution only in the range s±- < 
s(t) < s\ + is drawn in (a). The analytical result is in good agreement with that of the numerical 
simulation. 
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FIG. 4: Phase diagram on the (tJ 1 , Hq) plane. Transition line is independent of tJ 1 and is given 
by H = H c . 



SRM and SBM on the (tJ , Hq) plane is given by Hq = H c . The phase diagram is given in 
Fig. |U Furthermore, the long time average of s(t), (s(t)), depends discontinuously on Hq 
at Hq = H c as shown in Fig. El These behaviors are quite different from those of the DPT 
case driven by periodically oscillating field, F(t) = h cos(Qt) . The transition point h c 
for a fixed f2 between SRO and SBO depends on the frequency fl, and (s(t)) is a continuous 
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FIG. 5: (s) vs Hq for ry = 5. One finds that (s) is discontinuous at Hq = H c 



function of h. 



III. MFPT THROUGH THE CHANNELS 



We hereafter discuss the dynamics for H Q slightly above H c . Let us first consider the 
behavior obeying the equations 



s = s — s + eH , (e = + or • 



(14) 



for Hq > H c , i.e., F(t) is fixed to be either +H or —Hq. Equation for H Q > H c is 
integrated to yield 



In 



(s - s f 



si + s e s + s 2 -l 



2(3 S 2-1) s 2 + SeS + s 2_ 1 (S -S e ) 2 



6s F 



2(3^ -i) y /3g=Z 



arctan 



2s + s £ 



— arctan 



2s + s e 



(15) 



where s = s(0) and s e has been defined in Eq. (JUJ). Figure El displays three orbits given by 
Eq. ()15j) with sq = S- and e = +, which shows that s(t) approaches s + in the limit t — > oo. 
One observes that s(t) stays for a long time in the vicinity of s = — l/v3 for Hq slightly 
above H c . The small region including the position s = —l/y/3 is called the 'channel'. From 
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FIG. 6: Three orbits of the equation of motion (|14j) with e = + for Hq > H c . The values of Hq 
are set to be Hq = 0.386 (dotted line), 0.4 (dashed line), and 0.5 (solid line), where all the initial 
conditions are chosen as so = s_. 

the symmetry of the system, there also exists the channel near s = l/y/3 for F(t) = —Hq, 
as shown in Fig. E^b). Let us express the positions s c h of the channels as 



Sch 



-1/V3, ifF(t) = +H 
+1/V3, if F(t) = -H 



(16) 



The characteristic time is then defined as the time span that the state point s(t) passes 
through one of the channels for a constant F(t), either +H or —Hq. r c h can be estimated 
by integrating Eq. (jl4j) around s ~ s c h as follows. First, consider the case F(t) = —Hq. By 
setting u(t) — s(t) — s c h and assuming \u\ <^ s c h, Eq. (JTlj) is approximated as 



u 



-3s ch u 2 - (H - H c ). 



(17) 



This can be integrated to give 

u(t) = 



Hp — He 

3s ch 



tan 



y/3sch{H - H c )t 



with the initial condition u(0) = 0. r c h is estimated by the condition u(r c h) 

C n 



Teh 



(H -H c y/i- 



c 



(18) 

oo and thus 
(19) 
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FIG. 7: Time series of s(t) for (a) H = 0.388 and (b) 0.385 by keeping 77 = 200. 

Let us next consider the process that the state point s(t) passes through the channels 
under DMN. Figured shows temporal evolutions of s(t) numerically obtained for Ho = 0.388 



and 0.385. One finds that the time of passing through channe 



s increases as H approaches 

H c . The MFPT r through channels was calculated in Refs. [14, 3| by analyzing the master 
equation. In the present section, we will derive MFPT in terms of the time scales Tf and r c h 
from a phenomenological viewpoint without use of the analysis made in Refs. jlJ . Q] . 

The condition for passing through a channel is that F(t) continues to take the identical 
value either +Hq or — Hq for time longer than t c h- For H$ satisfying 77 > Tch, we obtain 
p{ T ch) — e~ Tch ' T f ~ 1, which implies that F(t) almost always satisfies the condition for 
passing through the channel. Therefore, f in the case of 77 > r c h is nearly equivalent to r c h, 
i.e., 

f ~ _ C ___^ . (20) 



(H -H c y/f 

In the case of Tf <^ r c h, on the other hand, Eq. (J3| gives p(r c h) 1. This fact implies 
that the probability that F(t) continues to take the identical value for time longer than r c h is 
quite small and hence that f is much longer than r c h because it needs a long time to satisfy 
the condition for the state point to pass through the channel, f in the case of r/ r c h is 
explicitly determined as follows. For a long f, let us divide f into subintervals each of which 
has the time span r/. The divided individual time series are approximately independent of 
each other. Therefore, 17/f is the probability that the state point passes through a channel 
once because f is MFPT through the channel. On the other hand, p{j c h) is identical to 
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FIG. 8: The dependence of f on Hq — H c determined numerically for rt = 5. f for a given Hq 
is obtained by integrating Eq. (jSJ until a time T and simultaneously by counting the number N 
of times passing through the position s = 0. Then f is evaluated as f = T/N. One finds that 
lnf -1 linearly depends on (i^o ~~ He)" 1 / 2 in the region of < ~~ H c -C 1, where the solid line is 
expressed as lny = Ax + i? with fitting coefficients A and B. 



the probability for s(i) to pass through the channel once by definition of the probability. 
Therefore we get the relation p(r c h) — t//t, which leads to 

C 



t 1 ~ Tj 1 e 



Ch/Tf 



Tj exp 



'T f (H -H e )W 



(21) 



with the constant C defined in Eq. f!19jb This expression agrees with the result obtained in 



Refs. 



14, 



151 ] . Equation (|21|) reveals that MFPT through the channel depends on H — H c 



in a stretched exponential form for tj r c h, and is quite different from the asymptotic form 
([20)1 . The above dependence of f on H — H c is confirmed in Fig. |HJ 



IV. PHENOMENOLOGICAL ANALYSIS 



In order to discuss statistical characteristics of the dynamics passing through the channels 
for Tf r c h, we here develop a phenomenological approach. The behaviors of s(t) for which 
we attempt to model are first summarized. The initial condition of s(t) is set to be in 
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FIG. 9: Time series of F(t) schematically indicating the time that s(t) jumps from s+ to s_. s(t) 
jumps at the time t* in this case, where F(t) takes the value — Hq longer than r c h for the first time 
in the time series of F(t). 

the vicinity of s + . If a time interval of F(t) satisfying the condition F(t) = —Hq becomes 
longer than r c h for the first time, then s(t) passes through s c h and approaches s_ in the 
time interval. See Fig. The event in which s(t) jumps from s + to s_ occurs only in this 
case. It should be noted that the jumps from s(t) > (s(t) < 0) to s(t) < (s(t) > 0) are 
approximately independent of subsequent jumps. 

Let us discretize the time t in the form t = kAt, (k = 1, 2, 3, ■ ■ •) as a simple approach 
to develop the phenomenological analysis according to the process noted above, where At 
is a certain small time step. Then, r c /j is discretized as r c ^ = n c hAt with the corresponding 
integer n c h- F(t) is assumed to keep the same value for the interval At, which is denoted as 
Fk = F(kAt). The conditional probability p that Fj+i takes the same value as Fj is given 
by 

p = e -At/r /j (22) 
and the probability q that F^+i is different from is therefore given by 

q=l-p. (23) 
The system is analyzed phenomenologically as follows: 

• We introduce the variable Sk at a discretized time kAt which takes two values ±1. 

• Sk and F k are initially set to s = +1 and F Q = +H , respectively. 

• Sk jumps from +1 (—1) to —1 (+1) only if F^ continues to take the identical value 
— H Q (+Hq) for a time interval longer than n c hAt. 
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Sk does not jump from +1 (—1) to —1 (+1) even though Fk continues to take +Hq 
(—Ho) for any time interval longer than n c hAt. 



A. MFPT r through the channel 

We first derive the exact expression for the MFPT f through the channel with the phe- 
nomeno logical approach. In considering the time series having Fk, f is evaluated as 



0<fc<Z 



(24) 

q=l-p 



Here g^i^ * s the number of the time sequences {Fj} for < j < I satisfying that Fj changed 
its value k times in each {Fj} and Sj jumps from +1 to —1 for the first time at t = lAt. 
Equation (J24)) is, furthermore, rewritten as 

f=TQ nch (q,p) (25) 

q=l-p 

with the differential operator T and the quantity Q nch (q,p) defined by 

d_ d_ 

' dq ^ dp 



T = At[q—+p—) (26) 



and 

Qn ch (^P)^J2 E 9 { kT ] <l k P l - k - (27) 

/ 0<fc<Z 

One should note that the q- and p-dependences in Q Uch are crucial and that q and p are 
considered to be independent in Eq. (|27)l . 

The explicit form of Qn ch {q, p) is then determined so as to satisfy the following conditions: 

• In considering any length of time series giving Fk, there exists a time interval of length 
n c h in the last of the time series, where all the Fk take the same value — H , i.e., the 
condition that Sk jumps from +1 to —1 is satisfied. 

• The condition for Sk to jump from +1 to —1 is not satisfied before the last time 
interval. 

One should note that the equality Q n (l — p, p) = 1 holds for any n, because the time interval 
described above always exists somewhere in a long time series. Particularly, for n = n c h, 
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Q nch {l — p,p) is obviously equal to the probability that Sj changes its sign, which must be 
unity for H > H c . 

As shown in Appendix |XJ the explicit form of Q n (<?,p) is given by 

(1 — p)qp n ~ l 



Qn(q,p) 



(28) 



[1-p) 2 -q 2 (l-p n - l Y 
where the condition Q n {l — p,p) — 1 is easily confined. Applying the operator (|2*U|) to the 
explicit form (}28|) with n = n c h yields the relation 

2 — p n cfc— i 



T = TQ nch (q,p) 
= At 



q=l-p (1 — p)p n ch-^ 

2 — e At / T f e - T ch/Tf 



(1 _ e-At/r^gAt/r/g-Wr/' ( 29 ) 

where the last equality is obtained by using Eqs. (J2*2~j) and (|2~3"j) with the relation r ch = n ch At. 
The exact expression of f is finally given by 

f = T f (2e Wr ' - 1) (30) 

in the limit of At — > by keeping r c h constant. Equation ()30|) qualitatively agrees for 
W 7 / > 1 witl1 the result (1211 . 



B. Distribution function -P(r) for the passage time r 

The distribution function P(r) for the passage time r through the channel s c h is deter- 
mined by solving the equation 



P(r) = 5(r-f)Q nch (q,p) 



9=1 -P 



(31) 



where £(x) is the delta function. The Laplace transform C[P](z) should be calculated in 
order to solve Eq. (j3Tf . By using the series expansion of Qn ch (q,p) given by Eq. (|2*7jl . the 
Laplace transform of P(t) is obtained as 

/•CO 

£[P](z) = / e-« P(r)dr = e- zf Q nch (q,p) 



q=l-p 



(32) 



9=1-P 



= E E ^V^V^) 1 "* 

Z 0<fc<Z 

— — y j_ ^ 

Equation (J32)l implies that C[P](z) can be obtained by replacing q and p by e~ zAt q and 
e" zA V in Qn ch {q,p), respectively, i.e., 



£[P(t)](z) = Q n Je- zAt q,e-^p)\ 



q=l-p 



(33) 
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Substituting the explicit form (|28|l for n = n c h into Eq. (|3*3*)l yields the equation 

( TtZ + 1 ) e -(^+rf 1 )T oh 

C[P(r)](z) = lT/Z + lje --x— (34) 

in the limit of At — > by keeping constant. By applying the inverse Laplace transform 
to Eq. (|34p. the distribution function P(r) is analytically evaluated in the series expansion 
as 



P(r) = r/e-^^^WO^-^cosh^ 



fc=0 



x=(t fc+1 ) 2 



. . , . . „. , to sinh(t 2 ) „, x t-i 2 coshf^) — U sinhftq 
9(h) cosh(ti) - fl(t 2 ) + 



t 4 3 sinh(t 4 ) — 3t 4 2 cosh(t 4 ) + 3t 4 sinh(t 4 ) 
^ + 



2 '• " 8 

(35) 



where ifc(r) = (r — kr^/rf and 0(t) is the Heaviside function defined by 

f 1 for t > 
= { ~ • (36) 

I for t < 

For details of the derivation of Eq. ([35]). see Appendix IB1 

Let us suppose to truncate the expansion ()35)1 at k = k c for an arbitrary fc c . It should 
be noted that Eq. (j3*5Jl gives the exact distribution for < t < k c r ch even though the 
truncation is executed, since all the terms individually include and so the terms for 

k > k c do not contribute to P(t) for r < A; c r c /j. The analytical result ()35|) is compared with 
the numerically evaluated distribution in Fig. One observes that the phenomenological 
analysis quantitatively explains the statistical property of passing through the channels. The 
characteristics obtained from the figure are summarized as follows. 

• There exists a region where P(t) =0 for r < r c h, which presents the minimal time of 
passing through the channels. 

• P(t) decreases exponentially for r 3> r c h, P(t) oc e~ aT with a constant a. 

• The rate a increases as Tf is increased. This tendency is consistent with the fact that 
the probability of passing through channels increases as r/ is increased since DMN will 
often continue to take an identical value longer than r^. 



15 




11)00 20OO 3000 r 4000 




2000 4000 6000 z 8000 



FIG. 10: Distribution function P(t) of time passing through channels. The theoretical results 
(solid lines) are compared with the numerically evaluated distributions (symbols) for (a) ry = 130 
and (b) Tf = 50. The strength of DMN is set to be Hq = 0.3852 in both of the numerical 
simulations. r c h corresponding to the applied DMN is estimated to be T c h — 134 via the condition 
-P(t) ~ for r < r c ^, where both (a) and (b) give the identical value of r c /j. The expansion in 
Eq. (|35|) is summed over for k < 60, i.e., k c = 60. 

On the other hand, the expansion disagrees with the correct value in an exponential 
way for r > k c T c h- Let us try to obtain the asymptotic solution of P{j) for r ^> r c h. 
Equation (jB2j) is approximated as 

1 



-z(r-T ch ) 



) 



for \z\ < t,} 



ch i 



(37) 



1 + (r - T ch )z 

where f is MFPT given in Eq. ()30|) . The inverse Laplace transform of Eq. (|3T|) is straight 
forwardly calculated to give 

1 IT- T ch 



Pt) 



T - T ch 



exp 



T - T ch 



(38) 
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which reveals that P(t) decreases exponentially with the damping rate a = (r — r c ^) 1 for 

T > Teh- 

C. Fourier spectrum determined by the phenomenological analysis 

We derive the Fourier spectrum of a time series s(t) by the phenomenological analysis 
to focus on the dynamical characteristics in the SRM phase. The Fourier spectrum I x (u>) is 
defined by 

^^h^M'i^df V (39) 

i.e., the ensemble average of the Fourier transform of a time series x(t). 

Let us first consider s (t) = sgn[s(£)]. Then the time series s (t) is expressed as 

s (t) = (-I)"" 1 , fort n _!<t<t n (40) 

with n > 1, where t n denotes the nth time to cross zero for s(t). Hereafter, t is set to be zero 
without loss of generality. By identifying that r„ = t n — t n _i is independently distributed 
according to Eq. (|35p. one obtains the Fourier spectrum of so(t) by the phenomenological 
analysis shown in Appendix in the form 



T (n) 4 y ^-^A 4 )| 2 U]] 

So{ ' fu 2 \1 + (e-^) y ™ 2 |1 + ( e -^\i2 



where 3?(X) represents the real part of X, and hin^^oo ■# = (Tn) = t is used. 

Substituting the explicit form of (e~ luJTn ) given in Eq. (}3"4*|) with 2 = iuo into Eq. ()41|) 
yields 

/4tA cj 3 t) + (4 - e- 2T ^ T f)ur f - 2e' T ^/ T f (u Tf cos ujt cIi + 2 sinujT ch ) 
S0 ^~\tTJ) (4 + uj 2 t 2 ){uj 2 t 2 - 2uT f e- T ^ T f smur ch + e~ 2T ^ T f) ' ^ ^ 

The above result is confirmed by comparing with the numerically evaluated Fourier spectrum 
for the normalized time series So(t) in Fig. 1111 

Let us finally modify the phenomenological analysis which is compatible with the numer- 
ically evaluated spectrum of the original time series s(t) without normalization. Instead of 
Eq. (jUJ), let us define 

s(t) = (-l) n - 1 [l-a(t-t n - 1 )) for t n -t<t<t n (43) 
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FIG. 11: Fourier spectrum of so(t) obtained theoretically (solid line) and numerically (dashed 
line). Tf = 130 is set for both numerical and theoretical results, and the strength of DMN is set to 
be Hq = 0.3852 in the numerical simulation. r c /j = 134 estimated in Fig. ^] is used for Eq. (|42|). 

with n > 1, where a (At) incorporates the wave form of the time series passing through the 
channel and is assumed to be a (At) = for At > r c h- Note that by setting a(At) = 
also for At < r c h the result of original phenomenological analysis is recovered. As shown in 
appendix O the Fourier spectrum I§{uj) for s(t) as a modification to i« (&>) is obtained in 
the form 

1 + \cl(uj)\ 2 + 2£[a(w)] 



I s (u) = I so (uj)- 



where 



a(uj\ 



Tch 

n,j I a(t)e~ iut dt. 
'o 



By approximating as a(At) = 1 + |sch| for < At < r c h, Eq. (jUJ) reduces to 



^oH I — ^ + 



cos^r c/l 



(44) 



(45) 



(46) 



Equation (J46j) is confirmed by comparing with the numerically evaluated Fourier spectrum 
in Fig. E21 



18 




10 



10 10 W oj 1 



FIG. 12: Theoretically obtained Fourier spectrum (|46|) (solid line) with \s c h\ = l/\/3 is compared 
with the numerical result (dashed line). Tf = 130 is set for both numerical and theoretical results, 
and the strength of DMN is set to be Hq = 0.3852 in the numerical simulation. r c h = 134 is used 
for Eq. ()42j) which was estimated in Fig. I1UI 

V. CONCLUDING REMARKS 



In this paper, we used the model equation (jHJ) under the dichotomous Markov noise 
(DMN) F(t) with a finite correlation time in order to investigate the dynamics of the mag- 
netization s(t) of the ferromagnet system driven by the magnetic field applied in one direc- 
tion, where its strength is constant and only the direction temporally changes. It was found 
that the dynamics of s(t) show two kinds of motion, i.e., the symmetry-restoring motion 
(SRM) and the symmetry-breaking motion (SBM), which are respectively observed when 
H is above and below the critical value H c . The transition line between SRM and SBM 
was determined only by the strength of the applied DMN and is independent of the corre- 
lation time Tf of DMN. By observing the distribution functions of s(t) for SRM and SBM, 
the ensemble average of s(t) discontinuously changes at H = H c . These results are quite 
different from those in the system driven by a periodically oscillating field []| ■ 

We then discussed the mean first passage time (MFPT) slightly above H c and found that 
it depends on H Q — H c and r/ as 

C 



t — Tf exp 



r f {H -H e )W 
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(47) 



This anomalous characteristics has a form similar to that of the average duration between 
neighboring phase slips of the phase difference in the phase synchronization observed in 
coupled chaotic systems j^J. Furthermore, a phenomenological approach was proposed to 
analytically discuss the statistical characteristics for H > H c . By obtaining the probability 
Qn(Q,p) of s(t) passing through the channel, the MFPT f, the distribution function P(r) 
of the passage time r and the Fourier spectrum of the time series were obtained 

by using the phenomenological analysis. The statistics obtained by the phenomenological 
analysis were found to be not only in qualitative but also in quantitative agreement with 
the numerically obtained results. 

In closing the paper, it is worth noting that the effect of DMN on nonlinear dynamical 
systems is generally quite different from that of Gaussian noise and DMN produces a new 
dynamical response of the systems. It is highly desired to examine the statistical char- 
acteristics obtained in this paper in laboratory experiments, and also in other numerical 
simulations, e.g., Monte Carlo simulations. 
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APPENDIX A: EXPLICIT FORM OF Q n (q,p) 

First of all, notice that each sample path of {k = 0, 1, 2, • • • ) corresponds to a symbol 
sequence of {+,—}. Let w be a symbol sequence of {+,—}, which is referred to as a string. 
For given sets A and B of strings, let AB be the set of all strings expressed as ab for a £ A 
and b G B, where ab denotes the concatenated string of a followed by b. In the following, 
the set composed of only one string w will be simply expressed as w. Furthermore, let {A) n 
be the set of all strings expressed as 

axa 2 ---a k (Al) 
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with di G A and < k < n, where k = means the zero length string, including the case 
di = cij for i 7^ j. Then, with this notation, every possible sequence starting with + and 
terminating with successive — 's of length n (n > 2), which firstly appears in that string, 
can be expressed as 

n symbols 

S n = ( + ( + )oo- (-)n- 2 )oo+ ( + )oo C ^ - ^. (A2) 

For a given set S of strings, let us define its "probability" as a function of p and q by 

P(q,p;S) = Y,P*M (A3) 

where P w (p,q) = p k ( w ^q l ( w ) denotes the probability for a string w, and k(w) and l(w) de- 
note the numbers of pairs of identical symbols ( , ++) and different symbols ( — h, H — ) 

appearing in w, respectively. For given strings w\ and W2, obviously, the identity 

(q,p) = P Wl a(q,p)P W2 (q,p) (A4) 

holds with a being the first symbol in w 2 , and we call w\o and W2 a decomposition of the 
string w = WiW 2 - By noting that 

Q n (q,p) = P(q,p;S n ) (A5) 
and considering decompositions of each element in S n , Q n (q,p) can be expressed as 

Q n (q,p) = P(q,p;S' n )P(q,p;R)p n -\ (A6) 

where 

S' n = (+(+) 00 - (-) n _ 2 ) 00 + (A7) 

and 

R = +(+)oo - . (A8) 
Since each element of S' n can be uniquely decomposed into a multiple of elements in 

Si = +(+)oo - Rn-2+ (A9) 

and inversely every multiple of elements in S" uniquely corresponds to an element in S' n as 
its decomposition, we obtain 

oo 

P(Q,r,s' n ) = ^[p(?,p;W (Aio) 

3=0 

1 

" l-P(q,p;S»)- ( H) 
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Finally, P(q,p] R) and P(q,p\ are calculated as 



CXJ 

P{q,p; R) = gYV = -r^— 
£ r~i 1 — p 



and 



3=0 



n-2 



j=0 i=0 

which yield Eq. (J2E|) together with flSEJ and (jATll . 



1 — p 



n-l 



pY 



(A12) 



(A13) 



APPENDIX B: DERIVATION OF THE PROBABILITY DISTRIBUTION 
FUNCTION P(r) 



Let us evaluate the inverse Laplace transform of Eq. (|34j) in order to derive the explicit 
form of P(t). For simplicity, the time is rescaled as t/ = 1. Then Eq. ()34|) is rewritten as 

(1 + z )e~( 1+z ^ h 



(1 + z) 2 - (1 - e -( 1+ ^)' 



which reads 

£[P(r + r ch )] = ( e -*(^-^)) = 

1 + 2 

Repeatedly applying the formula 



l-e _ ( 1+z ' T c;j 
1+1 



~ (l_ e -(i+^) n 

fl + z )2n+l ■ 



n=0 



c- 1 



-{l + z)T ch 



with /(r) = £ one obtains 

-(l+2)r ch ^ 



c- 1 



(1-e- 



(1 + ^) 2n+1 



-2n 



(2n)! 



where 0(r) denotes the Heaviside function Eq. (jSEJ), and the formula 



-i 



1 + * 



T 



m— 1 



(m - 1)! 



e- r ^(r) 



(Bl) 



(B2) 



/0) = e~ Tch f(r - r ch ) = e - Tch e~ Tch ^ f(r) = e~ Tch e~ Tch ^ X~ l [f{z)\ (B3) 



(B4) 



(B5) 



for any positive integer m was used. Thus, the inverse Laplace transform of Eq. (|B2|) reads 
P(r + r ch ) = e~ T <* }2 y 1 ~ e^e-^) e - T 9(r) 



n=0 



-(r+T, 



^1 - — -6{t), 



n=0 



(2n)! 



(B6) 
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where the identity 



e Tch e Tch dr q T = g r g Tch d T 



was applied. By noting the identity 



1 — g Tch dT 



y~v (— e Tch dT) k d k x r 

k=0 



k\ dx k 



E 

A;=0 



A'! 



x 



x=l 



Eq. ()B6|) is further simplified as 



oo oo 



P(t + r« 



c/i J 



n=0 fe=0 

feT „i (-1)* / d 



e -(r+r ch ) e -kT ch7I 

n=0 

oo 

(r+r ch ) ST^ 



-D* / d \ k r 2n 



k\ \dx J (2n)! 
d \ ^ 

— ) coshfv^xVW 
ax / 



fc=0 

.°° . r-r 2 ^ / d \ k 

{T+T Kh )S^„-kT Kh 4:ar--\\ T ) I a \ 



x=l 



x=l 



k=0 



k\ \dx J 



— I cosh \fx 



e -(r+r ch ) 9(t fe )L^-^ C osh V > 



k=0 



k\ dx k 



x=ti 



(B7) 



(B8) 



(B9) 



where tk(r) = r — kr c h and the formula ^^1 T^yy = cos h was used. After the replace 



ment r — > r — r c /j and the rescaling of time asr^r/rj and r c /j — > r c /j / r/ , Eq. ()B9j) reduces 
to Eq. lj5S|) . 



APPENDIX C: DERIVATION OF THE FOURIER SPECTRUM 

In this appendix, we derive the Fourier spectrum of the time series of the magnetization 
s(t) for the phenomenological analysis. Let Ti, t 2 , t 3 , • • • be a sequence of mutually indepen- 
dent random variables having the probability density P(r) obeying the condition P(r) = 
for r < r ch , i.e., r k > r ch . We introduce the variable 

s(t) = (-l) n [l - a{t - t n _0] for t n _! < t < t n , (CI) 

where t„ = Yl k =i T k an d a,(At) is a function satisfying a(At) = for At > r c ^. We assume 
that the time series of the magnetization s(t) is approximately expressed by s(t) with an 
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appropriate form of a(t). The Fourier transform of s(t) follows 

S(t)e- iu)t dt = ^(-l)"- 1 / [1 - a(t - t n ^)\e- iut dt 

n=l J 

N 

n=l 

TV 

= £(-i) n 



i — 1 — iut n — i 



-iujt 71 



i~n rr ch 

e~ iuJt dt- / a(t)e~ iuJt dt 
Jo 



a[uj\ 



n=l 



where 



a(u) = l- iw / a{t)e~ lu}t dt. 
Jo 



(C2) 



(C3) 



Considering the absolute square of Eq. (|C2j) and then taking the ensemble average, we obtain 



tN 



s(t)e~ lu>t dt 



2 N 



5>-*"»-a( w )| 



n=l 



+23R 



J] (-l)"- m e -^ E ^+i Tfc (e~ iu)Tn - a(u)) (1 - e~ iu)Tm a*(uj)) 

\<m<n<N 



(C4) 



and 



iV 



5{t)e-' luJt dt 



1 + \a(uj)\ 2 - 2^[(e- tuJT -)a*(u)] 



+2iV _1 3? 



^ ( e -^)«—i (( e -^) - a(w)) (1 - (e^ rfc )a*(w)) 



_l<m<n<N 

Since |(e _laJTn )| < 1, in the limit of iV — > oo, the last term in Eq. ()C5|) reads 



(C5) 



-23? 



J](-l) fe (e- iWr ") fc {{e- iulT ") - a{uj)) (l - (e'^a*^)) 



k=0 



-23? 



(e~^ T ") (1 + |a(cu)| 2 ) - a(o;) - (e-^ T ") 2 a*(o;) 
1 + (e~ iuJTn ) 



(C6) 



which leads to 



lim — 



s(t)e- luJt dt 



1 - (e~ iWT ") 
1 + (e~ iaJT ") 



(l + |aH| 2 + 23?[aH]). (C7) 



Noting that limjy^oo = (r n ) = r, we obtain 



1 + |a(^)| 2 + 29?[a(u;)] 
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where 



h (u) = 



OJ 2 T 



{e- iUTn ) 



1 + ( e - iu)T ") 

corresponds to the case that a(At) = for all At > and thus s (t) 



(C9) 

'-l) n for t n _! < i < 
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